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Abstract: In this paper, the theoretical expressions of wood thermal conductivity in the choral and radical di¬ 
rection are derived from the micro-structure of wood by applying some basic principles in physical mechanics. 
The thermal conductivities of about twenty species of trees were calculated by means of the expressions and 
compared with its experimental values under the same condition. The average relative error is about 5%, so the 
calculation result is satisfactory. 
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Introduction 

Thermal conductivity is required in scientific re¬ 
search and heat treatment of wood. Because wood is 
a kind of natural polymer organism whose chemical 
element and microstructure are very complicated, it’s 
very difficult to study its thermal conductivity theo¬ 
retically. Therefore, most scholars at home and 
abroad study presently the thermal conductivity by 
means of the experiments (Cheng et al. 1985; Koll- 
mann et al. 1968; Gao et al. 1993). Various empirical 
equations calculating the thermal conductivity of 
wood have been put forward based on the experi¬ 
mental data. But the question is that the empiric 
equations of the wood thermal conductivity derived 
from experimental data according to mathematics 
fitting have different functions, different adaptable 
range (species trees, moisture content, temperature 
and so on), with a larger errors, and lacking the 
theoretical foundation. In order to solve the 
above-mentioned problems, the universal theoretical 
expressions of wood thermal conductivity are derived 
in this paper in terms of basic principles of physical 
mechanics and wood physics. 

Mathematics deduction of thermal conduc¬ 
tivity of wood 

The carrier of energy in heat conduction of the 
general insulation solid is the transversal wave and 
longitudinal wave, i.e. the acoustic wave, caused by 
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the elastic vibration of the crystal lattice. The lattice 
vibration shows wave-particle duality. 

In order to derive the theoretical expressions of 
thermal conductivity of wood, we may assume that 
the phonon is the carrier of energy in heat conduc¬ 
tion of wood. The phonon shows the particle property 
of lattice wave, and moves irregularly in all directions 
in wood at sound velocity of u as gas molecules do, 
colliding each other and exchanging their energy. A 
phonon with velocity u, total energy E and mean 
free path L can transport energy E when it 
passes through the distance of L. 

Considering the double direction movement of 
various phonons, we find that energy is transported 
from the region of larger E to the one of smaller E , 
that is, from the region with higher temperature to the 
region with lower temperature. 

Now let’s consider the number of the phonons 
which penetrating vertically through one plane unit 
area S located in x=x 0 , which is shown in Fig. 1. 



Fig.l The derivation of thermal conductivity 

Supposing that the velocity u of the phonon re¬ 
mains unchanged and the number density of phonon 
per unit volume is n. The result of a statistical analy¬ 
sis shows that the number of phonons passed within 
the solid angle dQ(corresponding to the angle 0 and 
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0+d0 relative to the direction of ox axis, shown as 
a shaded ring in Fig.1) is in proportion to 
27r(rsin0)rd0 (i.e. the area of round ring AB). 

Therefore the number of phonons penetrating the 
unit area S per unit time in a small angle interval be¬ 
tween the angle 6 and 6+ 69 is 


n — ——— m cos d = — nu sin 6 cos ddd (1) 
Anr 2 2 


It is clear that the preceding collision of the pho¬ 
nons which are able to penetrating though the plane 
s from the left or the right side takes place in 
x = x 0 -Lcosd or x = x 0 + Lcos6 and then move 
freely though the plane S located in x=x 0 . We shall 
assume that the average energy of the phonon in x 
and Xo is E(x) and E{x 0 ) respectively. When the value 
of Lcos0 is smaller, E(x) may be expressed as fol¬ 
lowing 


capacity at constant volume. 

According to the above-mentioned mechanism of 
heat conduction of insulation solid, C v is just De¬ 
bye’ molar heat capacity, that is 

C v = 3RD(^L) (6) 


where D(-y-) is Debye’ function, 0 D is Debye’ 
temperature, that is 


0 = ^ 
D k 


f 


3 




(7) 


where h is Planck’ constant, k is Boltzmann’ con¬ 
stant, V is the volume of wood’s unit crystal cell, u 
is the sound velocity in wood. 

Substituting (6) into (5) we obtain 


d£ 

E(x) = E(x o ±£cos0) = £(x o )±£cos0 — (2) 

dx 


l =—Rup 0 LD\ 
M 


Or 


( 8 ) 


where the positive symbol corresponds to the pho¬ 
nons located in the right of the plane S and the nega¬ 
tive symbol corresponds to the phonon located in the 
left of the plane S. 

Multiplying formula (1) by formula (2) we can get 
the energy of the phonons which remove from the 
left (or the right) to the right (or the left) of the plane 
S along the direction between 0— 9+69 . The inte¬ 
gral of the difference in energy of the two ones be- 



transfers from the left to the right through the unit 
area S in per unit time 


equation (8) may be used to calculate the thermal 
conductivity in theory, but it is difficult to decide the 
mean free path of the phonon L. Generally L is 
related not only to the direction of heat current, but 
also to the frequency of sound wave in wood, tem¬ 
perature and density of wood, etc (Qian 1962). 
Therefore, it is very difficult to decide the expression 
of the mean free path L of the phonon theoretically. 
According to the experimental results and theoretical 
analysis, we assume that there exist the relationship 
between the mean free path L and temperature T , 
wood density, Debye’ characteristic temperature 
0 D as follows 


Q = -nuL — p cos 2 0 sin 969 = --nuL —— (3) 
6x 3 6T dr 


6T 6x 


According to the macroscopic law of heat conduc¬ 
tion, the quality of heat that transfer though the unit 
area in per unit time is 


Q = 



(4) 


After comparing equation (3) with equation (4), we 
obtain 


1 r 1 T Po r 

A - —utiL — = —uL—— C v 
3 dT 3 M v 


(5) 


where p 0 is the mass density of dry wood , M is the 
molar mass of wood cell, and C v is the molar heat 


£ = C, 



exp 

( I -\ 

C r A 

k 0 ° , 


l » P ° 


(9) 


where p A is density of air (values in normal state), Ci 
and C 2 are two constants to be decided (C* is of 
length’s dimension), and can be decided according 
to the experimental data of the thermal conductivity. 
On the basis of higher precision, we can use Ci 
=10.15X 10' 10 m, C2=7.987 (for heat current in the 
choral direction), and C, =12.60X10" 10 m, C^=5.283 
(for heat current in the radical direction). 

In order to simplify calculation, we can expand 
Debye’ function in expression (8) into Taylor’ series 
at 7=273.15 K, and expand T^ 2 in expression (9) 
into power series of t. Then, by substituting known 
quantity, we can obtain respectively simplified for- 
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mula of wood thermal conductivity in the choral and 
the radial direction 


\ ho 


= (9.257+0.0519) xl(T 3 



(W - m _l • K -1 ) 


( 10 ) 


l rad = (11.49 + 0.0645r) x 10 


& 


exp, 


' 6.0 X 


[Mj 


(W-m -1 -K -1 ) 


( 11 ) 


Comparison of theoretical values with ex¬ 
perimental values 

The relationship between thermal conductivity of 
wet wood (whose moisture content is w%) and that 
of dry wood is shown as the following equation 
(Kollmann etal. 1968) 

X w = (l + 1.25W%)X ( 12 ) 


where X may calculated by means of the equation 
(10) or (11). It is interesting to note that p 0 in equa¬ 
tions (10) and (11) is density of absolutely dry wood. 
In fact, the density p w with different moisture content 
can be turned into absolutely-dry density p 0 by the 
following equation (Cheng etal. 1985). 


1 + 0.01/TTV 

Po=-~ Pw 

1 + 0.01W 


(13) 


where K v is average volume shrinkage coefficient of 
wood. Evidently if the sound velocity, density, tem¬ 
perature, moisture content and volume shrinkage 
coefficient of wood are known, the wood thermal 
conductivity in the choral or the radial direction can 
be worked out. We have calculated the thermal con¬ 
ductivity in the choral of 23 species of trees at differ¬ 
ent temperature and with different moisture content, 
and compared it with experimental values under the 
same conditions, as shown in Table 1. 


Table 1. Theoretical and experimental values of thermal conductivity in chord direction of 23 species of wood 


Species of frees 

Moisture 

content 

W(%) 

Temp. 

ICC) 

Contract 

Coefficient 

k v 

Density 

pw 

(kg»m' a ) 

Sound velocity 
u(m»s‘’) 

Thermal conductivity 
(W-m' 1 K 1 ) 

Error 

(%) 

Theo. 

Expe. 

Pinus Koraiensis 

11.0 

16.5 

0.459 

456 

5057 

0.1080 

0.0989 

+9.2 

Cunninghamia lanceolata 

17.0 

23.5 

0.421 

459 

5079 

0.1166 

0.1070 

+9.0 

Pinus massoniana 

14.6 

24.2 

0,486 

490 

4958 

0.1235 

0.1221 

+1.1 

Podocarpus imbricatus 

13.6 

17.4 

0.419 

529 

4248 

0.1351 

0.1500 

-9.9 

Larix olgensis 

14.5 

13.9 

0.554 

702 

4658 

0.1627 

0.1605 

+1.4 

Pautownia fargesii 

10.5 

17.2 

0.334 

253 

4394 

0.0740 

0.0721 

+2.6 

P. fortunei 

13.0 

21.7 

0.320 

246 

4386 

0.0755 

0.1756 

-0.1 

P. catalpifolia 

12.7 

22.4 

0.344 

312 

4442 

0.0892 

0.1849 

+5.1 

P. tomentosa var. inlingensis 11.3 

17.6 

0.333 

321 

4360 

0.0889 

0.1896 

-0.8 

P. tomentosa 

13.0 

22.4 

0.261 

341 

4489 

0.0943 

0.0930 

+1.4 

Alniphyllum fortunei 

13.9 

24.5 

0.423 

450 

4609 

0.1186 

0.1210 

-2.0 

Tilia mandshurica 

10.9 

16.0 

0.447 

470 

5008 

0.1107 

0.0989 

+11.9 

Juglans mandshurica 

14.5 

16.1 

0.486 

481 

4973 

0.1166 

0.1105 

+5.5 

Catalpa duolouxii 

17.3 

26.0 

0.368 

482 

4688 

0.1271 

0.1233 

+3.1 

Melia azedarach 

15.3 

24.5 

0.413 

486 

4532 

0.1283 

0.1314 

-2.4 

Betula platyphylla 

13.2 

16.7 

0.494 

583 

4957 

0.1354 

0.1233 

+9.8 

Acer davidii 

12.3 

24.5 

0.388 

616 

4310 

0.1552 

0.1535 

+1.1 

Acer mono 

11.7 

25.0 

0.544 

659 

4471 

0.1631 

0.1686 

-3.3 

B. alnoides 

15.3 

24.2 

0.541 

674 

4879 

0.1627 

0.1698 

-4.2 

Fraxinus mandshurica 

12.1 

26.0 

0.548 

702 

4882 

0.1658 

0.1524 

+8.8 

Quercus mongoiica 

12.7 

25.0 

0.555 

721 

4288 

0.1809 

0.1710 

+5.8 

Homalium hainanense 

12.0 

24.0 

0.565 

900 

4494 

0.2107 

0.1965 

+7.2 

Quercus acutissima 

14.8 

21.0 

0.578 

963 

4433 

0.2265 

0.2396 

-5.5 


We have also computed the thermal conductivity in 
the radial direction of 18 species of trees, and com¬ 
pared it with experimental values, as shown in Table 
2 . 

It is clear from Table 1 (the thermal conductivity in 


the choral direction) that the maximum deviation from 
the experimental value is 11.9%, and the average 
error is 4.8%. Simultaneously, it is evident from Table 
2 (the thermal conductivity in the radial direction) that 
the maximum deviation from the experimental value 
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is 13.7%, and the average error is 5.3%. 

It will be seen therefore that the results of the 
comparison between the theoretical values and the 
experimental ones of thermal conductivity of various 
species of wood under different conditions is satis¬ 


factory. Therefore, it follows that the theoretical ex¬ 
pressions of thermal conductivity in the choral and 
the radial direction in this paper can basically reflect 
the objective law of thermal properties of wood. 


Table 2. Theoretical and experimental values of thermal conductivity in radial direction of 18 species of wood 


Species of trees 

Moisture 

content 

W(%) 

Temp. 

t('C) 

Contract 

Coefficient 

Kv 

Density p„ 
(kg»m') 

Sound velocity 
u(m»s~’) 

Thermal conductivity 
(W»m'' K' 1 ) 

Theo. Exper. 

Error 

(%) 

Pinus koraiensis 

11.7 

16.5 

0.459 

438 

5057 

0.1123 

0.1163 

-3.4 

Pinus massoniana 

15.0 

24.2 

0.486 

531 

4958 

0.1428 

0.1256 

+13.7 

Paulownia fargesii 

10.2 

20.6 

0.334 

250 

4394 

0.0757 

0.0791 

-4.3 

P. fortunei 

13.3 

21.7 

0.320 

262 

4386 

0.0805 

0.0884 

-8.9 

P. catalpifolia 

12.7 

22.4 

0,344 

327 

4442 

0.0961 

0.1000 

-3.9 

P. tomentosa var. tsinlingensis 

10.6 

17.6 

0.333 

306 

4360 

0.0884 

0.0907 

-2.5 

P tomentosa 

12.7 

22.4 

0.261 

335 

4489 

0.0967 

0.0942 

+2.6 

Alniphyllnus fortunei 

14.1 

24.5 

0.423 

418 

4609 

0.1192 

0.1314 

-9.3 

Tilia mandshurica 

10.6 

16.0 

0.447 

416 

5008 

0.1166 

0.1128 

+3.4 

Juglans mandshurica 

13.4 

16.1 

0.486 

481 

4973 

0.1243 

0.1279 

-2.8 

Catalpa duolouxii 

12.6 

26.0 

0.368 

472 

4688 

0.1300 

0.1233 

+5.4 

Melia azedarach 

15.1 

24.5 

0.413 

420 

4532 

0.1457 

0.1407 

+3.6 

Betula platyphylla 

12.6 

16.7 

0.494 

596 

4957 

0.1496 

0.1512 

-1.0 

Acer davidii 

11.5 

24.5 

0.388 

589 

4310 

0.1622 

0.1803 

-10.0 

Acer mono 

12.0 

25.0 

0.544 

670 

4471 

0.1821 

0.1942 

-6.2 

B. alnoides 

15.0 

24.2 

0.541 

752 

4879 

0.1850 

0.1954 

-5.3 

Fraxinus mandshurica 

12.1 

26.0 

0.548 

680 

4882 

0.1776 

0.1768 

+0.4 

Ouercus mongoiica 

11.5 

25.0 

0.555 

637 

4288 

0.1774 

0.1942 

-8.6 
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